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Abstract 

We use the results of integrable field theory to determine the universal amplitude 
ratios in the two-dimensional Ising model. In particular, the exact values of the 
ratios involving amplitudes computed at nonzero magnetic field are provided. 
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Universality is one of the most fascinating concepts of statistical mechanics [|I]]. Briefly 
stated, it says that physical systems with different microscopic structure but having in 
common some basic internal symmetry exhibit the same critical behaviour in the vicinity 
of a phase transition point. The point is best illustrated considering the singular behaviour 
of the various thermodynamic quantities nearby the critical point. For a magnetic system 
exhibiting a second order phase transition the usual notation is 
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where r = a (T — T c ) (a positive constant), h is the applied magnetic field and the limit 
towards the critical point r = 0, h = is understood. C, M, \ and £ denote the specific 
heat, the magnetisation, the susceptibility and the correlation length, respectively. 

The critical exponents a, /?,.. are the same for all systems within a given universality 
class and are related by the scaling and hyperscaling relations in d dimensions 

a = a' , 7 = 7', v = v' , 

7 = /3(<J-l), a = 2-2/3-7, 2-a = dv, 

a c = a/(35, j c = 1-1/8, v c = v/f38. 

The critical amplitudes A, B,.., on the other hand, depend on the scale factors used for r 
and h and are nonuniversal. However, universal ratios of amplitudes can be constructed in 
which any dependence on metric factors cancels out. Together with the critical exponents, 
these ratios further characterise the given universality class. The standard amplitude 
combinations considered in the literature are 

a/ A', r/r', (i) 

R C = AY/B\ = A 1/d £ , (2) 

R x = TD.B 6 - 1 , R A = A c D-^B- 2 /P , Q 2 = (T/T c )(U^r /u ■ (3) 
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A substantial progress was made over the last years in the derivation of nonperturba- 
tive theoretical results in two-dimensional statistical mechanics and quantum field theory. 
The solution of conformal field theories (CFTs) || |]] provided an almost complete clas- 
sification of universality classes for second order phase transitions in d = 2. In particular, 
it solved the problem of the exact determination of the critical exponents. The critical 
amplitudes, however, carry information about the scaling region outside the critical point 
and are not determined by CFT. In this respect, the Zamolodchikov's observation that 
specific perturbations of the critical point lead to integrable off-critical theories || is of 
crucial importance. The integrable theories obtained in this way, regarded as quantum 
field theories in 1 + 1 dimensions, are characterisable through the determination of their 
exact S-matrix. A number of physical quantities can then be computed using different 
techniques [0]. In particular, the results provided by the thermodynamyc Bethe ansatz 
(TBA) |7| and the form factor approach || |9|, [K| |ll[ enable the determination of am- 
plitude ratios and have been used for this purpose in the problems of self-avoiding walks 
PI and percolation |l3j (see also fi"4fl ). 



It is the purpose of this note to illustrate the derivation of universal amplitude ratios 
from (integrable) field theory through the very basic example of the two-dimensional Ising 
model. Of course, the purely "thermal" ratios ([I]) and (|2|) are exactly known for this case 
since the seventies, when the correlation functions of the Ising model at h = were first 
computed on the lattice |15], [T6f1 . The possibility to determine the ratios @, on the 



contrary, relies on the more recent realisation that the scaling limit of the Ising model at 
t = and h ^ is an integrable theory ||. 

We will regard the scaling limit of the two-dimensional Ising model as described by 
the euclidean field theory defined by the action 

S = Soft — t J d 2 xe{x) — h J d 2 xa(x) , (4) 

where Soft stays for the Ising critical point conformal action and a and e denote the 
magnetisation and energy operators, respectively. All the critical exponents are de- 
termined by the scaling dimensions x a = 1/8 and x E — 1 of the operators a and e: 
a = (d- 2x £ )/(d - x e ) = 0, (3 = x a /(d - x £ ) = 1/8, 7 = 7/4, S = 15, v = 1, a c = 0, 
7 C = 14/15, v c = 8/15. The physical dimensions of the two couplings in (f|) are r ~ m 2 ~ Xe 
and h ~ m 2 ~ Xa , m being a mass parameter. What is important for us is that the theory 
(£D becomes integrable when at least one of the two couplings is set to zero || . 

Choosing the scale factors for r and h amounts to fixing a normalisation for the 
conjugated operators e and a. We will proceed to determine the critical amplitudes 
within the standard CFT normalisation defined by the asymptotic conditions 

(a(x)a(O)) = \x\- 1/4 , x^O 
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(e(ar)e(O)) = |x|" 2 , x^O. (5) 

Let us begin with the amplitudes for the correlation length. The latter can be defined 
in different ways. For the time being, we will refer to the so called "true" correlation 
length defined through the large distance decay of the spin-spin correlation function 

e -kl/C 

(a(x)cr(0)) c ~ |x|(d _ 1)/2 , |z|->oo (6) 

where (• • -) c denotes the connected correlator. From the representation of the correlator 
as a spectral sum over n-particle intermediate states 

(a(x)a(0)y = 52\(0\a(0)\n)\ 2 e- E ^, (7) 

n 

it is clear that the "true" correlation length is the inverse of the total mass of the lightest 
state entering the decomposition above. For integrable models, the exact relation between 
the particle masses and the couplings appearing in the action is provided by the TBA 
[ |T^] . The result for the Ising model isQ 

m T = 2n \t\ , 
m h = C\h\ 8 / 15 , 

„ 4sinfr(l/5) /47T 2 r(3/4)r 2 (13/16)\ 4/15 

C = -. , x 5 ; ' '- , \, , ' = 4.40490857.. . 8 

r(2/3)r(8/i5) v r(i/4)r 2 (3/i6) J y 1 

Due to the invariance under spin reversal at h = 0, the magnetisation operator couples 
only to the states with odd (even) number of particles for r > (r < 0).When h ^ 0, a 
couples to any intermediate state. It follows 

£o = 2£, = 1/(2tt) , 

ic = l/C . (9) 

The specific heat diverges logarithmically in the Ising model (a = 0) and the specific 
heat amplitudes are accordingly defined through C ~ —A In r and analogous relations for 
A' and A c . Writing the partition function as Tr exp(— S) and the reduced free energy as 
/ = — 1/VlnZ, the specific heat per unit volume is given by 

C = = I d 2 x (e(x)e(0)) c ~ -2tt In mr , (10) 

or 1 J 



1 We denote by the subscript r (h) the quantities referring to the theory (||) with h = (r = 0). mh 
is the mass of the lightest among the 8 particles the mass spectrum of the r = theory consists of. 
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where (|5|) was used and ro is an ultraviolet cutoff. Recalling the relations between the 
mass scale m and the coupling constants r and h, the specific heat amplitudes in the CFT 
normalisation are simply 



A = A' = 2ir, A c = 2ir—. 

15 



The magnetisation per unit volume is given by 

Of 



M = — 



dh 



(11) 



(12) 



Vacuum expectation values in the CFT normalisation are exactly known for integrable 



models due to the TBA [17] and some more recent developments [18]. For the mag- 
netisation operator in the Ising model at r < 0, h = and at r = 0, h > they are, 
respectively^ 



(°)h 



2C 2 



15 (sin ^ + sin 



2- 



sm 



h i/i5 = 1.27758227.. /i 1/15 , 



15 / 



with 



B = 2 l/l2 (2n/e) 1/8 A 3/2 = 1.70852190.. 
we used the Glaisher's constant 

,1/2 



A = 2 7 / 36 vr- 1 / 6 



exp 



1 2 r 1 / 2 

- + - dxlnT(l + x) 
3 3 Jo V ; 



1.28242712.. 



(13) 



(14) 



(15) 



From (|13|) we deduce 



15 (sin ^ + sin \ + sin §) 
2C 2 



l 15 



0.0253610264.. 



The reduced susceptibility is defined as 

d 2 f 



X 



dh 2 



d 2 x(a(x)a(0)) c . 



(16) 



(17) 



The amplitudes T and V were computed exactly in |16[ integrating the h = spin-spin 
correlator. In the CFT normalisation they readf] 



T = 0.148001214.. , 
r' = 0.00392642280.. 



(18) 



2 In the following we will consider for convenience only positive values of h. Due to the simmetry about 
h = this involves no loss of generality. 

3 The relations one needs to pass from the lattice normalisation of to the one used here are 
a lat = 25/48 e i/8_4-3/2 a and [(T _ Tc)/Tc]lat = (7r/ ln (i + ^2)) T. 
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When t = 0, a is the operator which perturbs the critical point. The zeroth moment of 
the spin-spin correlator is then related to the vacuum expectation value as 

d*x (*(x)*{0))% = *' (a) h . (19) 
It follows 

T c = \d~ 1/& = 0.0851721517.. . (20) 

o 

The dependence on the operator normalisations (the only nonuniversal ingredient we 
had in our computation) cancels when the combinations (|l|), (0) and (|3]) are considered. 
The exact results we obtain for the ratios are collected in Table 1. The values for the com- 
binations of purely thermal amplitudes are well known. Concerning the ratios involving 
amplitudes computed at h ^ 0, to the best of our knowledge the only reliable extimates 
in d = 2 come from the series analysis of Ref. [19]. The value R x ~ 6.78 is quoted there 



together with a value for Q 2 which uses the "second moment" correlation length (see 
below). No previous result for Ra is known to us. 

We conclude this note computing the ratios which involve the correlation length am- 
plitudes using the second moment correlation length 



, 1 Jd 2 x\x\ 2 {a(x)a(0)Y 



^ 2nd ~ 2d f d?x (a(x)a(0)) c ' (21) 

The amplitudes (£o)2nd (£o)2nd could be exactly computed by numerical integration of the 
spin-spin correlation function at h = 0. Here, we will take a short cut with the purpose 
of illustrating a point which is particularly relevant for more general applications. In 
integrable models the spectral sum (0) is not a purely formal expression because the 
matrix elements it contains (known as form factors) can be computed exactly. While 
the Ising model at h = remains the single example for which the resummation of the 
spectral series is known, a remarkably fast convergence of the series emerged in the last 



years as a general feature of integrable models |10| , [12], |20|. This important circumstance 
makes the form factor approach extremely effective for obtaining accurate quantitative 
results at a relatively little cost. 

We compute (£o)2n<i and (Co)2nd in the leading form factor approximation (1-particle 
contribution for r > 0, 2-particle contribution for r < 0). Usingf] [^l[ 



(0\a(0)\9) T = (a) T , 

(0|a(0)|^ 2 ) T = i(a) T tanh^-^, (22) 



'Oil-shell momenta are parameterised as = (m T cosh (9, m T sinh( 
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one finds 



(£ ) 2 nd ^ 1/(27T) , 

(£o)2n d ^ l/(2Vl07r) . (23) 

(£c)2nd could be similarly extimated using the form factors for the r = model computed 
in ||20||. However, its exact value is immediately determined reminding that a is the 



perturbing operator at r = 0. Then the c-theorem sum rule holds 

c = ^[2ixh{2 - x a )} 2 J d 2 x \x\ 2 (a(x)a(0)) c h , (24) 

where the central charge c equals 1/2 for the Ising universality class. Combining with 
C P~9|) one gets 

{Zchnd = j-^ D l' m = 0.21045990.. . (25) 

Table 2 contains the results one obtains for the amplitudes ratios. The value (£o/£o)2nd = 
3.16.. is quoted in Ref. p3| as the result of the integration of the exact correlator; (£^2)2^ = 



2.88 ± 0.02 is the series result obtained in Ref. 19 
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Table Caption 



Table 1 . Exact amplitude ratios for the two-dimensional Ising model. The results for 
£o/£o' anc ^ ^ 2 re ^ er ^° ^ e "^ rue " correlation length. 

Table 2 . Results referring to the "second moment" correlation length with (£o)2nd and 
(Co)2nd computed in the leading form factor approximation. 
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A/A' = 1 
r/r = 37.6936520.. 

W£o = 2 
ifc = 0.318569391.. 

i?+ = l/v 7 ^ 
i? x = 6.77828502.. 
i? A = 0.0250658794.. 
Q 2 = 3.23513834.. 

Table 1 



(W£o)2nd - 3.162 
(i?+) 2nd ~ 0.3989 
(Q 2 ) 2 nd ^ 2.833 



Table 2 
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